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Q^ ■ Perturbative expansion in the nonperturbative confining QCD back- 

. ground is formulated. The properly renormalized as{R) is shown to 

Qh! be finite at large distances, with the string tension playing the role 

pi ' of an infrared regulator. The infra-red renormalons are shown to be 

^ ■ absent in the new perturbation series. 



^ : 1 Introduction 



The perturbation theory (PTH) expansion in QCD is an extremely useful 
tool, which has provided numerous links between theory and experiment [1]. 

The physical foundation of PTH in QCD is the celebrated asymptotic 
freedom phenomenon [2], ensuring the smallness of the physical coupling 
constant ag at large momenta (small distances). 

A formal extrapolation of the lowest order (one-loop) expressions for the 
renormalized to large distances (small momenta) is known to be plagued 
by the ghost pole, which signals an uncontrolled growth of the strength of 
interaction at large distances. It is believed that PTH expansion breaks 
down even before the Landau ghost pole is reached, so that the latter has no 
physical meaning, similarly to the case of QED [3,4]. 
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The impossibility of treating large distances in PTH of QCD strongly 
limits applications of QCD and from the formal point of view signifies the 
incompleteness of the theory. 

The logarithmic growth of the renormalized ag at large distances leads 
also to a basic difficulty, which is known as the infrared (IR) renormalon 
problem [5]. The IR rcnormalons appear as singularities on the positive real 
axis of the Borel parameter and destroy an attempt to sum up the PHT 
series. 

Thus the consistency of the PHT approach is seriously in doubt, when 
one expands around the free unperturbed vacuum. 

Alternatively it was understood during the last decades, that nonpertur- 
bative (NP) contributions are very important at large distances. In the frame- 
work of the QCD sum rules [6] the NP contributions haven been parametrized 
as local field vacuum averages using the operator product expansion (OPE) 
[7]. 

The phenomenologically found characteristics of the QCD vacuum clearly 
demonstrate that the NP fields in the vacuum are large (as compared to 
Aqcd). E.g. the gluonic condensate < (600Mey)^ > K^qcd [6]; 

lattice calculations yield even larger values [8] . 

The NP effects coded as power corrections in OPE [6] crucially change 
the behavior of the OPE series (as compared to the PTH series) when one 
is leaving the asymptotic freedom region and approaching large distances. 
To quote the seminal paper [6]: "The basic idea behind all applications is 
that it is the power terms (not higher orders in the series) that limit 
asymptotic freedom, if one tries to extend the short-distance approach to 
larger distances". 

It is specific to the local OPE, that power terms explode at large distances, 
demonstrating the importance of the NP effects but also the difficulty of the 
method to penetrate the region beyond the confinement radius. 

In this paper we propose to make a step forward in our strategy and 
instead of expansion in local condensates to consider from the beginning the 
PTH in the NP background of a general form. The framework for this exists 
for a long time [9-11]. Gluonic field is split into the background 5^ and 
perturbative gluon field and the new PTH is roughly an expansion in 
powers of g a^. The propagators of U/j, are exact Green's functions Gf^u in 
the background field 5^. 

A new element, crucial for the method, is that we are using a closed form 
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for G^u based on the so-called Feynman-Schwinger representation (FSR) [12]. 
This allows to display explicit dependence on B^, which enters through the 
parallel transporters P exp ig Jy B^dz^. In the limit of large Nc any term of 
the new perturbative series contains -B^ through the Wilson loops W(C, B). 
Averaging over an ensemble of background fields {B^} then reduces the graph 
to a product of averaged Wilson loops if N^. oo, and all the nonpcrturbative 
series is unambigiously defined by the expression for the averaged Wilson loop 
< W{C, B) >. At large distances the latter obeys the area law [13], and one 
obtains therefore a one-parameter description of NP dynamics in the new 
perturbation series. 

As a result one obtains a gauge-invariant PTH, even though the back- 
ground gauge fixing term is used in the Lagrangian together with Faddeev- 
Popov ghosts [14]. It is not a purpose of this paper to formulate Feynman 
rules for diagrams of the new PTH and to study systematically the renormal- 
ization procedure in the case of the NP background. Instead we concentrate 
here on the lowest one-loop contribution to the charge renormalization. In 
this way we attempt to answer the question: how the asymptotic freedom 
(AF) behaviour of cts, which is generated by the lowest gluon and ghost 
corrections in the background field, is maintained and how background pa- 
rameters enter into the final expression for a^. 

To do that one first must define the charge renormalization in a gauge- 
invariant and NP nonlocal formulation. We propose 3 definitions and two of 
them study quantitatively. 

We find that the NP background fields change the renormahzed charge 
in a drastic way for small momenta (large distances). 

The usual logarithmic growth of as{R) at large distances in the one- loop 
approximation disappears and as{R) saturates instead at rather moderate 
values. The mass scale m (or distance m"^) which governs this change of 
asymptotics is large, m ~ IGev, and physically is connected to the hybrid 
mass - or in the string terms - to the excitation mass of the transverse string 
vibration, ~ 27ro' ~ IGeV"^. 

As a result the asymptotic freedom can be kept at any distances in the 
Euclidean space and expresses itself as a saturation of as at as{max) < 0.5. 

We always use Eucledian metrics and do not analytically continue per- 
turbative expansions into the Minkowskian region, where physical thresholds 
appears. This will be a topic of separate publication. 

This property drastically changes the problem of the summation of the 
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PTH series. We show that the IR rcnormalons do not appear in the back- 
ground field, and the corresponding set of graphs has a finite sum. 

Hence the only obstacle left is instantons (instanton - antiinstanton pairs) 
, which require an explicit treatment . 

We stress that it not the purpose of our method to suggest an explicit 
way of splitting the fields into the background and pcrturbativc ones. Instead 
we choose from the beginning the NP background fields, they are considered 
as given and this input is parametrized, e.g. by the full set of correlators 
< F{xi)F{x2)...F{xn) >. 

The problem one has to solve in this case is how to find the perturbative 
series in this background, and we attempt to do it in the paper. Happily 
it appears that at large distances the only characteristics of the background 
which enters the answer is the string tension, and we take it to have the 
known phenomenological value. 

The plan of the paper is as follows In section 2 we introduce background 
field Lagrangian, gauge-fixing and ghost terms following the usual prescrip- 
tions [9-11] and write explicitly the gluon propagator. In section 3 three 
definitions of the running coupling are compared, and in Section 4 one of it 
- defined through the energy of static quarks, is studied in detail. Another 
definition, yielding the running coupling as{p) in the momentum space, is 
discussed in Section 5. 

Section 6 is devoted to the renormalons and the problem of summation 
of the perturbative series. 

In the conclusion we summarize results of the paper anl outline possible 
future developments. 

Three appendices devoted to the explicit derivations of formulas used or 
discussed in the text. In Appendix 1 we consider the contribution of the 
term Li, which directly connects NP background and the perturbative field 
a^, to the propagator of A^. 

Appendix 2 contains some detailes of derivation of the renormalization 
as{R) using the Wilson-loop approach. 

Appendix 3 contains a detailed disscussion of the one-loop graph in pres- 
ence of the confining background; most of the results used in the main text 
are based on the formulas explicitly derived there. 

In Appendix 4 we treat the gluon spin part of the gluon propagator in the 
background field. Finnaly, Appendix 5 is devoted to the explicit calculation 
of the IR renormalon contribution. 
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2 Background Lagrangian and propagators 



We split the gluonic field into the background and perturbative part 

A^ = B^ + a^ (1) 

As we stressed in the Introduction wc do not specify any principle of this 
division. For the purposes of the paper the exact properties of separation 
are immaterial, and we shall only assume that the NP fields 5^ ensure the 
area low of the Wilson loop < W{B) > with roughly the same string tension 
a as known phenomenologically, a ~ 0.2GeV. To dcvclopc the PTH in the 
background field B^ we use the well established background field formalism 
[9-11], specifically we extensively use method and notations of the paper [11], 
adjusted to the Euclidean space-time. 

The generating functional can be written as 

Z{J) Da det{—^)exp j d^x[L(a) - -{G^f + J>«] (2) 
where L{a) ^ Lq-\- Li (a) + ^2(0) + Lint{a) 

L2{a) = +^a^{IDl6i^^ - D^£)^ + igF^^)a^, = 

= lat[D'^Dl%,-D';^Df-gr''F^,]ai , (3) 



Lint = -^(i?^(S)a,-i?,(S)a^)«^r^^a>^-l(?V"'^<<r^^a^a^(4) 
The background gauge condition is 



G» = d,a; + gP'^Blal = {D.a.Y (5) 
and the ghost vertex [14] is obtained from ^ = {Dn{B)D^{B + a))ab to be 
Lghost = -9t{D,{B)D^{B + a))abeb (6) 
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The linear part part of the Lagrangian Li disappears if 5^ satisfies clas- 
sical equations of motion. We do not impose on this condition, and we 
show in the Appendix 1 that Li gives no contribution to the effect we are 
interested first of all - the modification of the AF logarithm. 

We now can identify the propagator of from the quadratic terms in 
Lagrangian .^2(0) — ^(G")^ 

Gf, = [Dl5,. - D,D. + igF,. + -^D.D.YJ (7) 

It will be convenient sometimes to choose C = 1 ^'^'^ ^nd up with the well- 
known form of propagator in - what one would call - the background Feyn- 
man gauge 

G''l = {bl-5^,-2igF,,)-' (8) 



3 Definitions of the running coupling 

In case of the free vacuum one can use different schemes of regularizations 
e.g. MS or AiS schemes ) and renormalize the coupling (7^, fixing the value 
of the vertex at some Euclidean momentum. For, example one may require 
that the quark (or gluon) propagator has the form of the free propagator at 
some momentum p, — —fx^ [15]. 

In the case of NP vacuum this kind of definition is not possible anymore 
and one should instead introduce a new NP notion of the regularized cou- 
pling. We consider here three examples of such definition . 

i) Definition of the running coupling through the interaction 
energy of heavy quarks (HQ) 

The interaction energy of a heavy quark and antiquark E[R) at a distance 
it! between them can be written as 

E{R) = E^p{R)-p^ (9) 

where Ej^fp is the purely NP part of the energy, while a.s{R) is a renormalized 
running coupling at distance R 

-.m = ^ = id + + - (10) 
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and go is the bare coupling, f{R) is to be calculated from the 1-loop diagrams 
in the background NP field. 

In case of free vacuum a^^^{R) has been computed to one loop in [16] 
through the bare coupling 

9'iR)=9o + 9'o^,C,lni^r (11) 

where 5 ~ 1/M is the UV cut-off. 

Renormalizing via MS scheme with the normalization mass /i one gets 

[17] 

"M5(/^, R) = «Ms{l + —[^{IwR + 7e) + T?;bo - ^Ca]} (12) 

77 Z iZ 6 

with bo = f iVc - , 7£ = 0.5777... 

The main feature of (11) and (12) is the logarithmic grouth of as{R) at 
large distances, and this feature produces both the Landau ghost pole [3,4] 
and the IR renormalon problem [5] . In the next sections we shall demonstrate 
that the strong confining NP background eliminates the logarithmic growth of 
as{R) and is this way solves both the ghost pole and IR renormalon problem. 

To calculate Ois{R) in the NP background we shall use the gauge- invariant 
definition 

E{R) = - hm din < Wc{B + a) >Z} (13) 

where contour C of the Wilson loop < Wc > is a rectangular R x T shown 
in Fig. 1. 

The averaging over fields S^, and regularization is assumed in (13) and 
will be specified in the next section. Note that E(R) (13) contains both PTH 
and NP contributions as in (9). We again stress at this point, that in our 
approach there is no problem of separation of a given E{R) into Enp and a 
PTH piece (which is of course not well defined) . Instead we assume that an 
ensemble of fields {-B^} is given to us, it is characterized by correlators 

{B,x,)...E,„.AB,X2)>B (14) 

and in particular by the string tension a, which can be expressed through 
the correlators. Our task is to construct the PTH in this given background 
and to express in particular as{R) through the nonperturbative parameters. 
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To perform renormalization in the background field formalism, one may 
use its gauge invariance [10]. One can visualize that all terms in (2) are gauge 
invariant with respect to gauge transformations "rotating" both and a^. 
It is convenient to split gauge transformations as follows [11] 

a^^^U+a^U, B^^B; + t\ (15) 

If we consider following [11] only amplitudes with no external lines 
(and we need only those), then renormalization of and 9 lines is inessential 
(corresponding Z factors compensate each other), while g, B,^ renormalize 
as 



^0 = Zgg , So = Z'J'B, Co = n (16) 
Now because of the explicit gauge invariance, Zg and Zb are connected [11] 

z;' = z'J' (17) 

as a consequence, the quantities gB^ and gF^^{B) are RG invariant, and 
therefore all characteristics of NP fields, like correlators (14) or string ten- 
tion arc RG invariant and should be considered on the same footing as the 
external momenta in the amplitudes. 

ii) Definition through the S-field two-point function 

This is equivalent to the method used in [11], where two 1-loop diagrams 
have been used of Fig. 2 and 3 with two external B lines. To make averaging 
over fields {B^j^} we multiply each diagram in the coordinate space with the 
parallel transporter in the adjoint representation, 

^{x,y]B) = Pexpig / B^dz^ (18) 

Jy 

and consider both gluon and ghost propagating in the field with propa- 
gators G^^{x,y, B) and G°-^{xy]B) respectively. The product 

H{x, y) ^< $^^(0;, y, B)r;^Gf,{x, y, B)G^\x, y, 5)1^ >b (19) 
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depicted in Fig. 4, yields the one-loop correction to the renormalized charge, 
which in absence of field reduces to the second term on the r. h. s. of 
(11). Using results of the next section one can again prove that dependence 
like ln\x — y\ disappears in (19) when |x — ^/| ^ m~^,m > IGeV . This will 
be clear from discussion in the next two sections. 

In a similar way one calculates two- and more loop corrections, using the 
same two-point diagrams as in [11]. 

iii) Definition through the finite box of size L. 

This definision is similar to that of refs. [18]. One fills the box with 
the NP field and considers one-loop correction to g^. As in [18] one gets 



f{L) = gl + ni{L)gt + ... (20) 



with ni — hQ — \hi, 



-§^ln detG ^In detG 

^0 = /ON ' hi = — - (21) 
9^i-o[B) a^^o{B) 

Here ?7 is a dimensionless parameter characterising strength of field B^, ro(-B) 
is the effective action (= | / F^^{B)d'^x). 

For the abehan (nonconfining) field used in [18] one obtains the asymp- 
totics [18] 

oo 

ni{L) -i^oo E(^n + sJnL)/L^ (22) 

n=0 

with r„, s„ some constants, sq = Thus one recovers the usual free- 

vacuum renormalization with the infra-red logarithmic growth. If instead 
one uses the confining field Bj^, the asymptotics of ni(L) changes and the 
logarithm InL disappears (the study of IndetG in case of confining back- 
ground was done in [19] for the case of nonzero temperature T, a similar 
analysis for T = will be published elsewhere). 

In the next sections we shall study in detail the first and second defini- 
tions of the running coupling and shall find its one-loop renormalization. As 
will be clear both definitions yield the same qualitative phenomenon: the 
running coupling as{R) saturates at large R at some finite value, and the 
mass parameter m, which cuts off the logarithmic growth of as{R), is of the 
same order in both definitions. 
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4 Calculation of as{R) through the Wilson 
loop formalism. 

In this section we shall follow the way outlined in the first definition of as{R), 
i.e. we shall use E(R) given by (13). To this end we need the PTH expansion 
of the Wilson loop < W{B + a) >, i.e. expansion of the powers in gluon op- 
erator a^. We start with the Wilson contour C of an arbitrary configuration 
and divide the contour into pieces Aa;^(n), writing identically 



W{B + a) = -^trPexpig J^{Bij, + a^^)dx^ 



AT 



^trP n [1 + ig{BMn)) + a^(x(n))) Ax^(n)] = (23) 



^ f N fc-1 

= TT*'^^ 1 11(1 + igB^Ax^) + n (1 + i9B^,Ax^)iga^{k)Ax^{k) x 

\.n=l n=l 

N 

X Yl {l + igB^Ax^{n)) 

n=k+l 

One can represent the sum (23) as a sum of Wilson loops iy(")(5) with 
n insertions of the field igafj,{k) on the contour C: 

oo 

WiB + a)^W{B) + J2{igrW^^\B;X{l)...X{n))dx,,il)...dx,M (24) 

n=l 

Now we must integrate over the fields Da^^ with the weight given in (2). In 
the lowest order < W{B) >a= W{B) for the properly normalized weight, 
and in the order 0{g^) one has a one-" gluon" exchange diagram shown in 
Fig. 5. 

-g^< W^'^ >a dxdy = -g' J $"^(x, y; B)GlZ'^^{x, y; B)^^'{y, x, B)dx^dy,{25) 

where are parallel transporters (18) with fundamental indices written 
explicitly, a, (3 = 1, ...Nc 

^sa,p, ^ t'X-^Gf, , (26) 
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and is the gluon propagator in the NP background field given in (7) 
and (8). 

It is appropriate at this point to specify the notion of the gluon in the 
background and to distinguish it from the free gluon, interacting only per- 
turbatively. We propose to call it the " hybrid-gluon" , or /i-gluon, since as 
we shall sec, in most cases it will aquirc a mass and behaves as a hybrid. 

In the next order, 0((7^), one obtains five diagrams shown in Figs. 6-10. 

Next step is ti average over the ensemble of NP fields {-B^}. It is not 
necessary to specify the ensemble and prescribe the measure; we can use the 
cluster expansion to obtain for the simplest term [13]. 

< W{B) >B= expY,^^ / da{l)...da{n) < F(l)...F(n) > (27) 

n=l 

Exact gauge invariant form of the correlator < F {!)... F{n) > is discussed 
in [13]. The lowest order term in (27) e.g. has two independent structures 
the Kronecker one and the full derivative term [13] 

/ < F^,{u)^u,v)Fp,iv)^v,u) >= i5^p5,,-5^J,p)D{u-v)+diDi) (28) 

It is D and the Kronecker structures of higher correlators that bring about 
the area law of the Wilson loop [13]: 

< W{B) >B= exp{-aSmin) , S large (29) 

where 

(7 = ^ y (fuD{\u\) + J < FFFF > +... (30) 

Thus any finite member of terms in the sum (27) or the whole sum, if it 
converges, yields the area law for the Wilson loop. In what follows we shall 
assume that the area law with a = (Tphen holds for all Wilson loops provided 
their size (both width and length) is large [13, 20], i.e. for the rectangular 
RxTin Fig. 1. 

R,T:^Tg (31) 

where Tg is the gluonic correlation length [20] , defined by the range of corre- 
lators D{u) etc. 

We stress at this point that the assumption (29) is very natural from the 
point of view of the cluster expansion [13], and moreover it is supported by 
the lattice calculations [21]. 
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With the help of (27) one easily obtains the first term in (9); using [22] 
one has 

rR roo rR roo 

Enp(R)^2R dX dvDi\v)+ XdX diy[-2D(X,iy) + Di(X,iy)] 
Jo Jo Jo Jo 

+higher correlators. (32) 

We note that at small R, R < Tg, [22,23] 

Enp(R) ~ CR"^ + O(R^) 
where 0(i?^) comes from higher order correlators, while at large R [22] 

Enp(R) ^aR-Co (34) 

We refer the reader to ref. [24] for detailed discussion. 

Now we turn to the /i-gluon exchange diagrams of Fig. 5-10. First we 
must regularize them, following the method of [25] . It was shown there that 
all divergencies of smooth Wilson loops are regularized though introduction 
of some minimal distance 6 on the contour C or by dimensional regulariza- 
tion and renormalized) by introduction of a common Z factor (heavy mass 
renormalization) and charge renormalization, generated by divergencies in 
diagrams of Fig. 6,7,9,10. 

We shall follow the same procedure as in [25], but we are interested not 
only in the divergent parts of diagrams as in [25], but also in the large R 
dependence with and without background field B^. 

We start with the simplest /i-gluon exchange diagram. Fig. 5, given by 
(25). The first step is to use for G/^i, the Feynman-Schwinger representation 
(FSR) [12] 

POD 

Gf^ix, y;B)^ dse^'iDzU^pix, y; B)f, (35) 

•/ 

where 

K=\l' zl(X)d{X); 
$^(x, y, B) = PpPexp ig f B^dz^exp 2 ig f F),p{z{X))dX (36) 

Jy Jo 

and differs from the usual parallel transporter (18) by the ordered (due 
to Pp) insertions of the gluon color magnetic moment interaction 2gFxp. 



(33) 
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The use of FSR is the central point of the method since it enables one to 
write explicitly dependence on 5^ - it enters via As a next step we 
must average (25) over the ensemble {B^}, i.e. to compute an average of an 
expression of the form. 

< ly(^) >s~< ^t^pt^ >B (37) 

where $j? is in the adjoint while $ in the fundamental representation. To 
understand better the structure of (37) let us take the limit N(. ^ oo and 
represent the adjoint hne in by the double fundamental hne, as suggested 
by 'tHooft [26]. We obtain two adjacent closed contours Ci and C2 as shown 
in Fig. 11. When both contours large, i.e. for R,T ^ Tg, and Nc —>■ 00 one 
has 

< Wc,{B) >B< Wc, >B +0(4^) ~ exp{-a{Si + S2)) (38) 

where Si, S2 are minmmal areas inside contours Ci, C2. 
Combining now (25), (35) and (38) we obtain 

-g''< W^^^ >B= -g^ / ds{Dz),ydx^dy,exp{-K - a{Si + S2))6^, (39) 

t/ 

(Note that we have omitted for simplicity the Fxp term in (36), which brings 
about spin-orbit corrections, see Appendix 4 for discussion.) 

Eq. (39) describes propagation of a scalar particle interacting with the 
contour via two strings, the world surface of which is given by 5*1, 5*2. 

We now show that for our choice of contour in Fig. 5,11, with the hori- 
zontal size T much larger than vertical size R, the resulting interaction tends 
to zero. Indeed for a given point ^(A) on the line of propagator, see Fig. 5, 
11 the interaction does not depend on the -24(A), Zi{\) which are in the plane 
of the contour Rx T, and depends only on the transverse coordinates 2:2, ^3. 
The interaction at the "time" A can be written as 



V = ai^zl + zl + zl + ^zl + zl + {T- z^y) - aT (40) 

where we have normalized to the case with no /i-gluon exchange. It is easy 
to see that for a generic case T — Z4, Z4 ^ Z2, Z3 (according to (13) we must 
finally take the limit T — > cx)) one has 

v^^^i±A + A±A)^o (41) 



13 



Thus we come back to the case of a free gluon exchange and obtain from 

(39) 



^ <w >B J J (x4 - 2/4)2 + ^2 4^ji y ) 

In a similar way we treat the diagrams with more gluon exchanges, like in 
Figs. 8, 12 and exponentiate the result (42), obtaining the lowest order term 
( — I ) in (9). One should mention that in multiple gluon exchanges like 
Fig. 12 there appear in general correlated gluon pairs, triplets etc. with 
average distance between gluons of the order of {a)~^/'^. These correlated 
exchanges are not taken into account during exponentiation and should be 
accounted for separately, yielding a Yukawa-type correction to the pertur- 
bative part of E{R) starting with the order 0((?^). We shall discuss these 
corrections in a separate publication, they do not change results of this paper 
concerning large distance logarithmic behaviour. 

Now we come to the fourth order diagrams, needed to renormalize 
shown in Figs. 6, 7, 9, 10. 

First we discuss the structure of divergencies and logarithms in case of 
no background B^. Referring the reader to the Appendix 2 for details, we 
can e.g. write the resulting integral for the diagram of Fig. 10: 



R Fig.io 167r^ J \x\\x — R\^ An R 6 

In dimentional regularization one would have 

R C 

In— — > h In/iR + const (44) 

d £ 

For Figs. 6 and 7 one obtains (again 5^ = 0) 

{^U,.e,r = (13 - ~ ^ Uir') - '-^C,ln^ (^5) 

where n('^^) is the trasverse self-energy part of gluon and ghost in the coor- 
dinate space: 

n('-=) - (46) 
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Integrating over d'^r we obtain the same logarithm as in (43). 
Some more care is needed for the graph of Fig. 9; leaving details to the 
Appendix 2 we mention only the final result 

f{R) _ jS-QC, R 

- 2nR ^^T ^^^^ 

Summing (43), (45) and (47) together we obtain 

fiR) = ^Inj (48) 

in agreement with (11) and (12). 

Consider now the case of nonzero i?^, and let us average each of the 
diagram W^^^{B) in Fig. 6, 7, 9, 10 over {B^}. We again use the limit Nc oo 
and replace each gluon line by a double fundamental line. Choosing the gauge 
fixing constant ^ = 3 we eliminate the logarithmic part (and divergence) of 
Fig. 9 and disregard it in what follows, since now we are only interested in 
the question: what happens with logarithms in (43), (45) and (47) when the 
NP confining field is included. 

Following the same line of reasoning as we used for the diagram of Fig. 
11 and imposing the confining field 5^ we obtain for diagrams of Fig. 6,7 
and that of Fig. 10 three surfaces in Figs. 13 and 14 respectively, filled 
with confining film due to the area law. As in (39), we must consider again 
horizontal size T of diagrams in Figs. 13,14 large and a result the surfaces 5*1 
and 5*2 do not produce any force on the gluon line when T — > oo. This is not 
true, however for , which provides a strong confining force between the 
fundamental lines which border 5*3. 

As a result we must consider the behaviour of the self-energy part Yl{x,y) 
in Fig. 13 , defined by the FS representation: 

ll{x,y)r^J^ dsj^ dsDzDzexp{-K - K - aS^) (49) 

It is easy to show ( for details see Appendix 3) that Y{{x,y) behaves as 
j^Zy)i small \x — y\, \x — y\ <^Tg, and has an exponential fall-off at large 

1^:^ - y\ 

JJ(a;, y) ~ exp{—mi\x — y|), \x — y\ ^ oo (50) 
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Introducing this cut-off behaviour in (46) one ol:)tains the following inte- 
gral which replaces the original logarithmic expression in (45) (see details in 
Appendix 2). 

i/n| m ^ f (i^(i? - r) + ^^(r - R)) (51) 

There are two regimes of behaviour oi I{R): the region of the asympthotic 
freedom : 

1 Re 

miR < 1, 1{R) ^ -In— (52) 
R 

and the infrared region 

rriiR > 1, 1(R) « - / — e""^!'' ^ -In (53) 

RJs r R mid 

The same situation occurs for the diagram of Fig. 14. Here the cut-off 
function exp^—mlR — x\) appears under the integral in (43), producing the 
same integral I{R) as in (51). 

As a final result of this section we formulate the new expression for f{R) 
taking the NP confining into account (modulo constant (nonlogarithmic) 
terms) 

^fNp{R) = (54) 

with I{R) defined in (51). 



5 Renormalization of as in the definition ii) 

In the second nonperturbative definition of ctg through the two-point func- 
tion H{x, y) (19), depicted in Fig. 4, we can again exploit the oo hmit. 

In this case each adjoint line, including ^{x,y, B), is represented as double 
fundamental line [26] and each closed fundamental contour Wc{B), after av- 
eraging over fields {B^}, acquires the area-law factor exp{—aSmin{C)) - in 
agreement with our assumption (29). Thus each of the diagrams of Fig. 4 
can be depicted as in Fig. 15, the difference between (a) and (b) is that for 
the gluon Green's function there are in addition gluon spin operators inserted 
on the gluon line, leading to the spin-dependent effects, which we disregard 
in the first approximation. 
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The diagrams of Fig. 4 can be written in the form similar to (49), which 
is discussed in detail in the Appendix 3 



o2 

H^,^{x, y) = j^bo{d^d^ - d'^Si,^)Y[{x, y) (55) 

where fl is defined as 

n(^,2/) = y^ dsj^ dse-''-''{Dz),y{Dz),ye-'' (56) 

with 

Y _ 0(1) I 0(2) I 0(3) CcyN 

and correspond to the minimal areas inside the 3 contours depicted in 
Fig.l5. 

We have omitted in (57) for simplicity the gluon spin terms in the Green's 
function G in (19). 

At small distances \x — y\ one can follow analysis of Appendix 3 and 
obtain that 

n(^^^)- (^_^)2.-4 (58) 

This leads to the famihar answer for the Fourier-transform of H^i, [11] 

2^ 2 
H^,t.{p) = -^^^ {Pt^Pv - P^^ixv) {Ne - In ^ + const) (59) 

with 

6o = ^A^c, A^. = --C + ln47r, 
3 e 

C is the Euler constant, C = 0.577. 

At large \ x — y \ one must take into account the confining dynamics 
present in V and find the corresponding mass. To this end one must look at 
the expressions (55), (56) more carefully. 

To understand the structure of H{x, y), one can dissect the graph in Fig. 
15 along the dash-dotted line and one obtains a system of two gluons 1 and 
2 connected by a pair of strings, as shown in Fig. 16. One of the strings 
goes through the fixed point 3 which appears due to the dissection of the 
straight-line parallel transporter $(a;, |/; B) - see Fig. 15. 
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One can show that the mass of the system of two gluons with two strings 
in the configuration of Fig. 16 is larger than the corresponding mass without 
the additional condition of string passing through the point 3: 

m(l,2;3) > m(l,2) = m2 (60) 

Thus we come to the conclusion that H describes the Green's function 
of a two-gluon glueball with a quantum nember J^'-' = 1 (and with an 
additional "heavy gluon" propagating along the straight line of ^{x, y, B), if 
we keep the point 3 in Fig. 16, and without it if we calculate a more light 
system according to the inequality (60)). 

An estimate of m(l, 2) with J^'-' — 1~+ has been done within the proper 
time Hamiltonian in [27] 

m{l,2)^2GeV ±0.15GeV (61) 

A similar estimate has been obtained in the lattice calculations, where usually 
two neighbouring members of the same multiplet L = 1,S = 1 are studied 
[28]. 

Using the arguments given in the Appendix 3, one can approximate 
Y[{x,y) as 

n(x.,).Q --^<-;i-^i' (62) 

which after renormalization yields the Fourier-transform 

n (p) ~ Co In !!i±^ + const (63) 

A comparison of (59) and (63) shows that up to power corrections two 
expressions have the same /i dependence, and hence the same Z factors. 
Actually all the change from (59) to (63) can be described as introduction 
of a new (fifth) external momentum component = m^^ so that the total 
momentum squared is in case of nonpcrturbativc background equal to p^+rn^. 
This fact does not change the Gell-Mann-Low equations which to one-loop 
order are [15] 

-fcoT^ (64) 
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and solving (64) with boundary conditions as in (55) and (63) (neglecting 
power terms as before for large (m^ +P^)a^o ^ 1) obtains 

= ^ (65) 

Introducing the A parameter (65) is rewritten as 

47r 

"-NP 

We have used notation Aatp to stress that this new A parameter may 
differ from the free vacuum h.QCD and needs a new determination from the 
experimental data. We stress here again that the appearence of m\ together 
with is possible because of renormalization invariance of gB^ [11], which 
in its turn yields invariance of a and m2 ~= 

At large p one obtains in (66) the familiar AF logarith. This logarithm is 
still large for p = if one chooses a low normalization parameter K^p^Kmp ~ 
^QCD- Thus the AF phenomenon can be maintained down to very low p in 
the Euclidean region. This is similar to the situation which we have observed 
in the Wilson-loop renormalization scheme, eg. (53). 

There, however, the mass mi corresponds to the one- Wilson-loop dynam- 
ics, while in our case actually a doubled film defines dynamics for m(l,2), 
with (Ta = 2(7, so that the mass m2 = \/2mi. This fact demonstrates that 
the IR cutoff is scheme dependent. (Note that our definition of the renor- 
malization scheme differs from the usual one, and the scheme dependence 
actually means the dependence on the physical process in question). 

6 Renormalons and asymptotics of the per- 
turbative series 

In this section we study the consequences of the modified AF (66) for the 
properties of the PTH series . We shall concentrate for definiteness on the 
Euchdean correlator of e.m. currents ]Y'"^'{Q'^) as in [29, 30] and shall write 



19 



it as in [30] in a form exhibiting the contribution of the renormalon [5] 



(67) 

where we have separated out of the PTH series the contribution of the set 
of diagrams shown in Fig. 17, denoted by Af]- We shall first discuss the 
infrared (IR) renormalon, [5], i.e. the contribution from the graphs of Fig. 17 
from the domain of integration over k with k <^ Q. If the coupling constant 
as is normalized at (we shall take it at + w? and (66) is used for a^, 
we have 

^^^-^^■^^^^Uo -Q^^^'^J^T^' ^^^^ 
The analysis of the integral (68) made in the Appendix 5, yields 



/ 

Jo 



C^'"<I7W) = | ,^ , „>2„„, ,<1 '^^^ 



with the notation no = In — 



As a result we obtain the following contribution of the graphs in Fig. 17 



where 

/ (n-1)! , n<2no 

We remind that the estimate (71) is an upper limit for q^. 

First we discuss a connection with the previous results [29,30] on the IR 
renormalons. In this case (no background field) and no — >■ cxo; g„ is 

defined as in (71) for any n and the Borel transform of AH is as in [29,30] 



An^U^^Vf^^^t^ (72) 

Stt 
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the Borel integral (72) exibits IR renormalon pole at t = |^ which makes the 
Borel summation of the PTH series prohibitive. 

In our case, when no is finite, the situation is drastically changed. Namely, 
having in mind that the estimate (71) is an upper limit of both for n > 2no 
and for n < 2nQ, we can take the sum in (70) which does not diverge and 
even docs not need Borelization. Leaving details for the Appendix 5 we quote 
only the final result 

A n = ^ E ^^^^^re-^^dt ^ '-^^ (73) 

There are no poles in the plane of the Borel variable, as can be easily seen 
from the absense of the factorial growth of the coefficient qn at arbitrarily 
large n. We thus conclude that the convergence of PTH series is radically 
improved when the strong nonperturbative background is taken into account, 
and in particular the IR renormalons disappear. 

We now turn to the ultraviolet (UV) renormalons [5] . Those come from 
the same scries of graphs, Fig. 17, but from the integration region with k"^ ^ 
Q^. It is easy to understand that the UV renormalons do not change due to 
the modified AF logarithm (66), since in this kinematical region k"^ ^ 
and the change is immaterial. Hence we have as in [30] the estimate 

^UV) = -±-^^a:i-^r-\n-l)\ (74) 

which brings a pole at negative value of the Borel parameter, t — —j^, and 
hence the Borel sum of the PTH series is well defined. 



7 Conclusion 

We have introduced in this paper a new perturbation theory, based on the 
exphcit separation of NP and perturbative fields, Eq.(l). 

The former are defined by their vacuum correlators as in (27), which at 
large distances enter only as one constant - the string tension a (30). 

This is the input of the theory. The PTH expansion in powers of pertur- 
bative field is then defined unambigiously in terms of the NP input and 
coupling constant ctg. 
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The main emphasis of the paper is the drastic change in the behaviour of 
renormahzed as at large distances as compared to the usual (free) PTH. 

Two different NP definitions of the running coupling constant are used, 
Eqs.(9,13) and (19), and in both cases the one-loop contribution is shown to 
be bounded at large distances , in contrast to the logarithmic growth of the 
free one-loop result. 

It is important that Gell-Mann-Low equations are not changed and the 
resulting behaviour of as{p) is characterized by only one new parameter, 
as in (65-66), which is RG invariant. The mass m depends on the process 
and in the typical situation of a gluon exchange in the e+e^ annihilation, is 
of the order of one GeV. As a result the magnitude of as{p) everywhere in the 
Euclidean region of is limited and the Landau ghost pole is not inevitable. 
Rough estimates show that as{p) is less than unity in the Euclidean region; 
enabling one to use perturbation expansion at all distances. 

As a consequence there appears a solution of the long standing problem 
of IR renormalons. It is shown in Section 6, that IR renormalons are absent 
when the running coupling has the form (66). Hence the Borel summation 
of the new perturbative series is not prohibited by the IR renormalons, and 
one encounters in the right Borel-parameter half-plane only the distant sin- 
gularities due to instanton-antiinstanton II pairs. 

We do not consider here important questions, e.g. how inclusion of II 
pairs into the NP vacuum fields influences the properties of the new pertur- 
bative series. 

We also do not discuss the behaviour of as{p) in the Minkowskian region, 

< 0, which shall be presented in a subsequent paper. 

During the earliest stage of this study the author intensely discussed 
with P. van Baal and G.'tHooft. Their useful remarks and suggestions are 
gratefully acknowledged. 
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A.Yu.Dubin, B.L.Ioffe, A.B.Kaidalov, L.N.Lipatov, M.I.Polikarpov and K.A.Ter- 
Martirosyan. 
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APPENDIX 1 



The contribution of the term Li 

Let us calculate the contribution of the term Li to the gluon propagator. 
If one denotes by <>„ the integral Da^, with the weight L{a) as in (2), we 
obtain 

j d%d'^y2Gp^^^{xi,yi)DpFp^^{yi)DxFxu2{y2)G,,^u,^{y2,X2) (Al.l) 

The gluon Green's function G^iy is given in (8) and depends on the back- 
ground field Bf^, as well as and Fp\. Averaging over {Bf^} is done in 
(A. 1.1) using cluster expansion . To get a simple estimate of A we replace 
G^u by the free Green's function G'^^} and take into account that 

< DpFp,MDxFx.,{y2) >b^ < Fp^MFx^M > (Al.2) 

^yip oy2\ 

and for the latter we use the representation [13] in terms of two independent 
Lorentz structures, D{yi — ^2) and Di{yi — y2). The contribution of D is 
(that of Di is of similar character) in the momentum space 

AS,,(fc) - (^''^MiM2-feMifeM2) ^(^) (A1.3) 

where 

D{k) = j d^ye'''yD{y) (A1.4) 

Insertion in (A.1.4)of the exponential fall-off for D{y) found in lattice calcu- 
lations [31] yields 

D(k) - < '''(Q) > . (Al 5) 

With II !^ IGeV, D{k) 0.12. Thus A^{k) is a soft correction fast 
decreasing with k as k~^ to the perturbative gluon propagator. 
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APPENDIX 2 



Calculation of the Wilson-loop charge renormalization 

Here we discuss in more detail calculation of graphs, corresponding to 
Figs. 6-10, which contribute to the charge renormalization in the "Wilson" 
renormalization scheme, Eq.(9). 

We start with the graph of Fig.6. A typical term contributing to the 
amplitude can be written as 



/(j (j 
A^i^(yi, x)—A^),{x - y)A^p{x - y)—Ax^,(y, ^2) x 

X dyiu,^dy2^j,^d^xd'^y (A2.1) 



Let the upper and lower horizontal lines in Fig.6 correspond to yu G 
[0,T],yii = and 1/24 ^ [0,T],|/2i = R, respectively. According to our 
discussion in Section 4, the first and the last gluon propagators in (A 2.1) 
are not affected by confinement when T — > 00 and one obtains in (A 2.1) 

j(6) ^ dx.^d. Hix, y%-S^. (A2.2) 

m \R-y\ 

Applying Fourier transformation to both propagators in (A 2.2) yields 

j{6) ^ y^MlW (A2.3) 
\R-f\ 

In case of no background field n(^) is simply a product of two free gluon 
propagators and one gets 

J"^'^ ^Lln^ (A2.4) 

where is the UV cut-off parameter. 

We introduce the following notation for the one-loop renormalization of 
as in the scheme (i) 

as{R) = af\l + f{R)af^ + ...) (A2.5) 
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In the case of no background the diverging contribution of graphs Fig. 6 and 
7 to f{R) is [25] 

(6,7) 13-3^ R 

/(i?)=^^^2ln| (A2.6) 

This should be compared to the momentum space calculation [17] of no 
background case (Feynman gauge, C = 1) 



(Q) = af\l + ^[(-^C,)(-- + C - ln47r + In \) + ^C,] (A2.7) 

4:71 6 e /i y 



a 



fr/ 5^,, 2 ^ , ^ 0\ 31 



where C = 0.577 in the Euler constant and in the MS scheme one should 
omit the factor (| + C — ln47r) in (A. 2. 7). We now consider the case of 
nonzero background field, and take n('") ^ estimated in Appendix 3, Eq. 
(A3.21) 



Insertion of Y{{r) (A2.8) into (A2.3) yields instead of (A2.4) 
where 

rOO rl'if roo rjj- 1 1 

^M = i Jd^^L -e-»"*(5e(«-r) + -e(r-«)) (A2.10) 

It is clear that the integral over dx converses at a; ~ 1 and therefore one can 
replace in the integral over r, rriichx — > m*,m* > mi. It is interesting to 
study two limiting cases: 1) m*R <^ 1 and 2)m*R 3> 1. One has 

1 Re 
I{R) -In — , m*i?< 1 
it d 

/(i?) f« , m*i?» 1 (A2.11) 

R m*d 

We turn now to the grphs of Fig. 8 and 9. Due to the nonabelian character 
of fields the graph of Fig. 9 is proportional to 

tr(tHW) = ( fN - — (iV^ - 1) (A2.12) 
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The first term on the r.h.s. of (A2.12) can be considered as the "abehan 
extrapolation'' of the graph of Fig. 8 to the crossed kinematics of Fig. 9, 
while the second term on the r.h.s. of (A2.12) represents the pure nonabehan 
contribution. In the case of no background one obtains therefore the following 
nonabelian contribution to the Wilson loop from the graph of Fig. 9 



f dyindy2ndxiudX2u { \o 



which gives the following divergent contribution to f{R) (cf(A2.5)) 

/^^'(i?) = ^^^lnf (A2.14) 

This should be compared with the momentum space calculation of [17] (we 
keep only divergent (logarithmic) term) 

Aa.(g) = N • (3 - In ^ (A2.15) 

To simplify our discussion we shall use now the gauge ^ = 3 and therefore 
one can disregard the graph of Fig. 9. This is true also in case of nonzero 
background, since we are interested here only in the modification of the 
divergent term In j due to the background. 

We now turn to the graph of Fig. 10. 

Its contribution to f{R) is 



R 167r2 'J \ x\\x-R, 

where | ^ enters from the static gluon propagator A(yi — x) while the 

favctor I x — ^ is from the triangle part r{Xjy2,y3) of Fig. 10. In case 
of the confining background one has to insert in the integral in (A2.16) the 
factor, exp{—mi \ a; — ^|), which yields again the function I{R) as in (A2.9). 
This factor appears in the asymptotics of T{x,y2,y3) due to the confining 
area law with the surface Ss, shown in Fig. 14. The asymptotics holds true 
when I a; — -R I is large and the mass mi is the same as in (A2.8), since at 
large \ x — y2 \ , \ x — y^ \ and fixed 1^2 — ^3 1 the function r(a;, 7/2, J/s) behaves 
as the two-point function Yl{x,y) discussed in Appendix 3. 
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Combining all contributions from all graphs of Fig. 6-10 one obtains 
finally 

(total) ^j^-^ = ^I{R) ,bo = ^N (A2.17) 

For small R , niiR <^ 1, one has due to (A2.11) the usual renormalization 
factor 

as{R) = + oi'-^^^^^T + ••■) (A2.18) 

ZTT 

while for large R, rriiR ^ 1, the function RI{R) tends to a constant and the 
unlimited growth of the renormalized charge as{R) is frozen. 
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APPENDIX 3 



Properties of the self-energy function Y[{x, u) 

In this appendix we study the properties of the self-energy part Y\{x,y), 
Fig. , both in the coordinate and in the momentum space. We consider 
for simphcity spinless particles in the loop, which is true for the ghost loop, 
inclusion of spin is discussed in Appendix 4. 

The FSR for n(^, v) is [12] 

roo POO 

Ui^:y) = J^ dsj^ ds{Dz),y{Dz),ye-''~''-^ (A3.1) 
where V reduces to aSmin for large loops [13] 



Introducing as in [24] " center-of-mass" and relative coordimates, and 

R^ = ^M>^) + ^^i^i^) ; r^ = z^{\)-z^{\) (A3.3) 
and new variables // instead of s, s 

2s 2s A* + A* 



Ui^^y) -T^r^. r $-.DR ■ Dr^^>'^^-^S:^^''^-^} (A3.5) 



we obtain 

djj, dji 
2^ Jo 2p? 

where T = \x — y\. The boundary conditions for trajectories Rn{T),r^{T) 
are: 

RM = y^; Rt.{T) = x^- r^(0) = r^{T) = (A3.6) 

as disscussed in [24], one can replace V — aS^in by the expression V — 
(T|r(r)|, with the accuracy of ~ 5% for the mass eigenvalues. In this case one 
can go over to the Hamiltonian form in (A 3.5) using the relation 

J DR Dr e- 1 ^'^'^ ^< R ^ x,r = 0\e- 1 ^'^^\R = y;r ^ > (A3. 7) 
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Since the motion in i?^ and r4 coordinates in free one can easily calculate its 
contribution and the result for Y[{x,y) is 



ajj,iajj,2e 2 



where G{rf,rin;T) is the Green's function of relative quark-antiquark (or 
gluon-gluon) motion, 

G(0, 0; T) = 5: C„(^)e-="('^)'^, (A3.9) 

n 

Cnifj) — |'0n(O)P,'0n ^ud £n(/i) are eigenvalues of the Hamiltonian 

//Vn = £nV'n; = |^ + (7|r1 . (A3.10) 

For linear potential Cn{jl) = |^ and does not depend on n , while is 

= (2/i)-V3^2/3^(^) (A3.11) 

where a(n) are numbers, obtained from the solution of the dimensioless equa- 
tion, with the asymptotics 

a{n)^{^-^f\n + lr/' , n ^ oc (A3.12) 

One may get a false impression that expressions (A 3.9- A 3.12) are nonrel- 
ativistic. In fact relativism is contained here in the p, variable which finally 
yields spectrum very close to that of the relativistic quark model [24] or 1 -|- 1 
QCD. 

Insertion of (A 3.9) into (A 3.8) yields 



where 



with 



l[{x,y) = AJ:J ^^e-^"('^-'^^)^ (A3. 13) 

M„(//i, /.2) = ^^^^ + Sniil) , (A3.14) 

A- ^ 
" 4(27r)7/2 
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At large T = |a: — |/| the asymptotics of (A 3.13) is obtained performing 
the diii,dii2 integrations via the steepest descent method [24], which yields 

Xli^^y)-^Y.^f^ne-^-^ (A3.15) 



where = 4-v/ct(^^)^/^ is obtained from the minimum of Mn{ni,ii2) in 
variables fii, fi2- 

Therefore at large T the correlator H decays exponentially 

expi-MpT) 
IK^'^) ^^375 

where Mq is the minimal mass value. 

At small T one should use (A 3.13) and take into account , that in this 
situation the sum can be replaced by the integral 



dne 



-A(A)T(n+i)2/3 



- l^lr^ ' ^ - (2/i)-^/V^/3(^)V3 . (A3.17) 
Insertion of (A 3.17) into (A 3.13) yields 

n(-.!/) = ip^f *2e-^ = jg^ (A3.18) 

This is the exact answer for the free scalar gluon loop diagram. It is 
remarkable that the sum over bound states reproduces the correct small 
distance behaviour. 

Finally we shall write the answer for the momentum space. Choosing 
p — (0, P4) one has 

At large n one can replace the sum by the integral of the type {D = 4+e). 

/ — = -- ^"("^0 - 'P^ ^) (^3.20) 

Jmo (M — la)^--^ £ 
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Integration over dO of the logarithmic term in (A 3.20) yields finally the 
effective logarithm which appears in (63) 



r'"^ cos^ ede ■ ln(mo - ipA sin 9)^^ ln(m^ + pi sin^ 9) (A3.21) 

J-Tv/2 4 



On another hand each term with a given n contribiites a simple pole in 
n(p), which can be separated out looking more closely at the integration 

2 

region near 6' ~ ±|. Indeed writing sin 6' = cos rj ~ 1 — ^ one can rewrite 
the integral in (A 3.19) as 



(5 + M„f)5/2 S^Mn' 3 

where S — — ipA- Finally the contribution of the poles is 



(A3.22) 



The sum (A 3.23) is diverging for large n, calling for the regularization 
which we have explicitly done in (A 3.20). 

The fact that n(p) has only simple poles equidistantly placed at —p"^ = 
= m^(ri + |) tells us that the difference Yiip) — Y[{Po) is proportional to 
the function, ip{z) = 

Indeed, one can write 

UiP) - UiPo) - ^(^-^) (A3.24) 

fit 

2 

which asymptotically at large ^ ^ 1 behaves as 

1 °° R 

V^(^) - In z - - V (A3.25) 

where i?„ are BernuUi numbers. In this way we again obtain the logarithmic 
term similar to (A 3.21). 

Finaly we derive one more form of Yl{x,y) which also was used in the 
main text. One can rewrite (A 3.15) separating in M„ the constant term, 
not depending on n: 

M„ = mo + M(n) (A3.26) 
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such that at n = Umin, M^rimin) = 0. 

Having in mind that dM^ = ni^dn, one can rewrite (A 3.15) as 

n(.. ,) ^ / .M^MV^e- . ^^M^fi:^ (^3.27) 

This is exactly the form used in the text, e.g. in Eq.(62), and in the Appendix 
2, Eq.(A. 2.8). 
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APPENDIX 4 



Gluon spin term in the gluon propagator 

The gluon propagator in the background field, Eq.(8) can be written in 
the form 

G={Dl.i + ^g±F)-' (A4.1) 
where 1 and are unit and spin matrices of the gluon respectively 

ig{Y.F)^.^-'2igF^u (A4.2) 

In the situation, when G enters in the gauge-invariant combination with two 
parallel transporters as in (26), (37), one can write the Wilson loop average 
as in (37) using the cluster expantion. E.g. for the adjoint Wilson loop 

< Wadj >=< ^f{x, y, ci)^F{y, X, C2) > (A4.3) 

where is defined in (36), and the paths Ci, C2 form the closed loop C = 
Ci + C2 one obtains [32] 

< Wadj >= expJ2^-^l A(l)...A(n) « F(l)...F(n) » (A4.4) 

n=l 

We have used the following compact notation 

(1) (2) 

A(A;) = da{uk) -J^^n + Yl ^fk (A4.5) 
with the condition, that e.g.the combination 

J^dn « F(l)...F(k)...F(n) :$>^J2drk « F{l)...F(z(rk))...F(n) » 

(A4.6) 

where z{rk) lies on the path Cj. At the same time the combination da{uk) is 
multiplied with F{uk) where lies on the surface Smin inside the contour C . 
Therefore keeping in A(A;) only the first term da{uk) one obtains the usual 
area law at large distances [13]: 

< ^adj >no spin— ^Xp{~'^adj Smin) (A4.7) 
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In contrast to that the terms with or generate perimeter-type 
contribution. Indeed the cumulants <^ F{l)...F{n) ^ are effectively nonzero 
when all coordinates are within the interval of the size Tg, vacuum correlation 
length, Tg « 0.2/m [31]. 

The operator X]'^*'' d^k enforces, as we discussed above, one of the coordi- 
nates of the cumulant to lie on the contour and therefore all others arc within 
the distance Tg from the contour, leading to the perimeter- type contribution 

< Wadj >spin^ exp{-(TadjSmin " hL) (A4.8) 

where b is some constant and L - the length of the contour C. 

For large contours C, which are studied in this paper, one can therefore 
neglect in tlie first approximation the gluon spin contribution as compared 
to the area-law term. 
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APPENDIX 5 



IR renormalons and background fields 

The integral in (68) can be rewritten as 

io + m2 io Jo ^ ^ 



where no = In — 



The sum in (68) can be now calculated easily (the interchange of the 
order of summing and integration is possible because the integral is in finite 
limits, and the sum is converging uniformly for i :< no). One has 

t{as{Q')t^r^'^^^^^^^ (A5.2) 



and the contribution of / in the sum over n is 
where 



bo .„2^ r° _e-^Hdt notdte-^ 1 



47r , Q^ + m^ 



since m > A. 

Insertion of (A 5.3) into (68) finally yields 



^Yl^^.L^bdm (A5.5) 



This result is used in Eq.(73) of the main text. 
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Figure captions 



Fig.l. The rectangular Wilson loop of the size R x T used to define the 
static quark energy E{R) 

Fig. 2. The gluon loop with external lines corresponding to the back- 
ground field B^, used in [11] to calculate P function. 

Fig. 3. The same as in Fig. 2 but for the ghost loop. 

Fig. 4. The nonperturbative generalization of the diagrams of Fig. 2 and 
3 with the full account of the background in the propagators of the gluon (a) 
and the ghost (b). 

Fig. 5. One-gluon-exchange in the Wilson loop between the points x and 
y. Indices a, S and (3, 7 are the color fundamental indices of the gluon at 
points X and y respectively. 

Fig. 6. The gluon-loop diagram for the gluon exchange inside Wilson 
loop. 

Fig. 7. The same for the ghost-loop diagram. 

Fig. 8. Two-gluon exchange diagram from iteration of one-gluon-exchange 
interaction. 

Fig. 9. Two-gluon exchange crossed diagram, containing both iteration 
of one-gluon exchange (abelian part) and the nonabelian part. 

Fig. 10. Gluon exchange with the three-gluon vertex in the Wilson loop. 

Fig.l 1. At large Nc the one-gluon-exchange diagram of Fig. 5 goes over 
into the product of two Wilson loops. 

Fig. 12. Multiple iterative gluon exchanges. 

Fig. 13. At large the diagrams of Fig. 6 and 7 contain three Wilson 
loops with areas 81,82,33. 

Fig. 14. At large Nc the diagram of Fig. 10 contains three Wilson loops 
with areas 6*1, 5*2, S3. 

Fig. 15. The diagrams of Fig. 4 contain at large three Wilson loops. 

Fig. 16. Cutting the diagram of Fig. 15 along the dash-dotted line reveals 
two gluons 1 and 2, connected by two fundamental strings. One of the strings 
passes necessarily through the point 3, lying on the parallel transporter $(-B), 
shown in Fig. 4 and Fig. 15. 

Fig. 17. Diagrams containing the renormalon contribution to the correla- 
tor of e.m. currents. 
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